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Abstract 

The bound-state spectrum of Hamiltonian H = H(\) is assumed real in a non-empty 
domain of multiindex A e j){physicai) _ rj^g name £ q Uan tum catastrophe (QC) is introduced 
for the scenario in which A leaves this domain in such a manner that an N— plet of the bound 
state energies degenerates at the boundary QD ( J )h y sical ) subsequently, complexifies. Our 
benchmark version of the QC scenario uses a family of crypto- Hermitian N by N matrices 
H (N \X). Up to the very QC horizon Qr£){p h y sical ) the explicit construction of the standard 
physical ad hoc Hilbert space = 'H^(X) (i.e., of a hermitizing metric = 9^(A)) 
is simplified and shown feasible at any N. 



1 Introduction 



The measurement of an N— plet of energy levels belongs among the most frequently encoun- 
tered forms of the experimental verification of theoretical conjectures concerning quantum 
systems. This is one of the key reasons why people often turn attention from a realistic 
Hamiltonian (where P stands for "physical" ) to its "friendlier" version yielding 
the same energies. This means that we can write, formally, 

H^^QH^Q- 1 . (1) 

Besides the most common cases in which f2 is required unitary, simplifications — > 
are accompanied by the loss of the manifest Hermiticity of the effective Hamiltonian, 
7^ At the first sight, this looks like a discouraging paradox. The non- 

numerical replacements of a complicated by its simpler avatar found just a 

few fully explicit realizations in the literature, therefore (cf., e.g., PQ). In parallel, the 
pragmatic numerical recipes of the form of Eq. ([1]) were more successful. In nuclear physics, 
for example, the use of non-unitary (often called Dyson) maps Q proved feasible and 
efficient [2J. 

Non-numerical theoretical activities in the field have been perceivably revitalized re- 
cently (cf., e.g., reviews [3j HJ 13 E|). People imagined that the friendly, non-Hermitian 
Hamiltonians with real spectra may be proposed and studied directly, without any 
recourse and explicit reference to their partners H^ p \ For this purpose, it appeared suf- 
ficient to translate, via Eq. (TjQ), the postulate of Hermiticity of = [H^y into the 
equivalent (sometimes called Dieudonn'e's [7] or crypto- Hermiticity [8]) requirement 

[HW] 1 Q = QH^ (2) 

involving merely the effective Hamiltonian. This relation also contains just the product of 
Dyson operators = Q^Q called, conveniently, Hilbert-space metric or, simply, metric [2j. 

Whenever the spectrum of energies remains real, the use of non-unitary maps Q does 
not lead to any conceptual difficulties (cf. Appendix A for more details). One simply 
has to keep in mind that the definition and/or an exhaustive description of a quantum 
system in question may employ either the manifestly self-adjoint Hamiltonian or its 
crypto-Hermitian reincarnation H^ F ' plus a sophisticated metric = ^ I. 
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In our present paper we intend to study and discuss what happens near the Kato's 
exceptional points [9], at which the spectrum ceases to be real [TO]. In particular, we 
intend to offer and describe a series of crypto-Hermitian quantitative models simulating 
the quantum catastrophes (QC) during which an N— plet of the bound state energies 
degenerates at an exceptional point and, subsequently, complexifies. 

2 The domains of hidden Hermiticity 

2.1 One-parametric Hamiltonians H(X) 

The difference between the use of the trivial metric 0^ = I and its more involved form 
q(S) -L J might look purely technical. During the first stages of the intensive development of 
the crypto-Hermitian representations of the operators of observables, this opinion prevailed. 
People felt addressed by the questions of non-emptiness of the domains of parameters in 
which the Hamiltonian was diagonalizable while its bound-state spectrum remained real, 

all E n (X) e R iff AG V {physical) . (3) 

The subtleties of the correspondence between Hermitian and simplified appeared 
inessential, in the weak-coupling perturbation regime at least (cf., e.g., [TTI [T2] ). 

A few years later, the emphasis has been shifted to the strong-coupling dynamical 
regime. People started paying attention to the existence and role of the Kato's exceptional 
points A = \ {EP) E &D^ hysical) [H E3j. Step by step, the studies of this possibility 
transcended the boundaries of quantum theory. The complexification of the spectrum, i.e., 
the often encountered fact that 

some E n (X) £ R when A i v {physical) (4) 

appeared relevant, e.g., in magnetohydrodynamics (where the related instabilities do exist 
and are measurable [14] ) or in the classical and experimental optics [T5] . 

The turn of attention to the models based on the use of nontrivial 0^ 7^ / changed 
the paradigm even inside quantum theory. In our present paper we intend to show that the 
new perspectives are truly opened by the possibility of explicit constructions of the metrics 
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in the strong-coupling regime. This challenging theoretical option is made important by its 
phenomenological appeal, connected either with the possible abrupt loss of correspondence 
PD at A — \( EP ^ or with the phase-transition- resembling violation (jlj) of the spectral reality 
beyond the EP horizon. 

2.2 Crossing the horizons Qj)(p h y aical ) 

In phenomenological applications of Eqs. ([1]) and fl2]) the quantum system is often studied 
in a fragile dynamical regime [16]. In the language of mathematics this means that the 
parameter A (typically, a coupling [U]) is located in a close vicinity of one of its EP values. 

Naturally, one has to distinguish between the parametric dependence of the complicated 
Hamiltonian = H^ p \\) and the combined parametric dependence of the simplified 
effective Hamiltonian = H^ F '(\) and of the effective metric = 0^(A). The main 
formal reason is that it is virtually impossible to speak about the operator H^ P '(X) in any 
nontrivial vicinity of \( EP \ Indeed, its Hermiticity is robust so that it cannot produce any 
non-real eigenvalues. In other words, this operator ceases to exist at any point A = A^ p ^ 
of the horizon. 

In contrast, there emerge virtually no relevant difficulties in the crypto-Hermitian rep- 
resentation of the systems near a quantum catastrophe. The dynamics using the doublet of 
operators \H^ F '(\), 0(A)] is flexible and can be prolonged till the horizon, i.e., up to the 
limit A — > \( EP \ In our recent papers [T7] - [25] on this subject, we solely paid attention 
to the study of the spectra in this limit. As a consequence, all of our results just concerned 
the geometry of the domain x*( phystcal \ In other words, we were interested in the behavior 
of the energies and not in the behavior of the wave functions and of their inner products. 
Thus, in spite of our initial optimism [20] , we only managed to parallel the Thorn's theory 
of classical catastrophes [26] by a rather formal combinatorial classification of the energy 
mergers [27] . 

In our present paper we intend to broaden our perspective and to fill the gaps. We 
will make use of the family of N by N matrix models as reviewed briefly in Appendix B. 
For the sake of definiteness, we shall merely select one of its most friendly one-parametric 
subfamilies. For these particular toy- model Hamiltonians H^ F \\), the construction of the 
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(up to now, missing) physical metrics 0(A) will be shown feasible and outlined in detail. 



2.3 The metrics 

The simulation of the catastrophe at A = must be based on our specification of the 

Hamiltonian and metric (cf. Appendix A). The usual, "friendly" inner product (f,gY F ^ = 
Yln=i fn 9n must be replaced by its "sophisticated" double-sum generalization 

N N 

(f,9) iS) = Y,Y,f™ e ™9n (5) 

m=l n=l 

defined in terms of a suitable physical Hilbert-space metric = > 0. On the level 
of abstract principles jl], the related amendment of the Hilbert space is 
of a vital theoretical importance since it changes the status of the Hamiltonian H(X) 
from "non-Hermitian" (i.e., unphysical) to "Hermitian" (i.e., observable). On the level of 
practical calculations, on the contrary, the interpretation of Eq. (J5J) as defining another 
Hilbert space THS S ' is redundant and may be, in the present context, treated as formal. We 
may always keep working inside where every use of the "correct" product (JSJ) may 
simply be translated, via Eq. (jSJ), into the friendly-space language in which the matrix 
©^ is inserted wherever necessary. This will allow us to keep using the Dirac's bra and 
ket symbols without any danger of misunderstanding. In particular, we may work with 
the usual Schrodinger equation 

H(X) |^-(A)} = Ej(X) l^-(A)) (6) 

in the Dirac's notation. In a useful additional convention of Ref. [I] the eigenvectors of the 
conjugate Hamiltonians may be denoted by the specific, doubled ket symbols, 

fft(A) |ifob(A)» = E k (X) |^(A))) (7) 

Such a convention simplifies the spectral representation of the metric (cf. Ref. [28] for 
more details), 

N 

©(A) = £>»(*)» «Vn(A)|- (8) 

71=1 

Once we now restrict our attention to the QC scenario in which an iV— plet of the bound 
state energies {E n (X)} stays real inside 2}( p/iys * ca/ ) an d complexifies beyond an exceptional 
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point A( £P ) G Qj)(s> h y sical ) ; we need not be particularly careful. Enhanced difficulties with 
the proper physical interpretation of the system should merely be expected due to the in- 
creasing sensitivity of both the strongly non-Hermitian matrix H = H(\) and increasingly 
ill-conditioned matrix = 0(A) to perturbations. Hence, in QC regime, a full-precision 
knowledge of metric 0^ would be highly desirable. This may be identified as one of our 
present main constructive targets. 



3 Solvable crypto-Hermitian toy models in the strong- 
coupling QC regime 

As long as iV < oo is finite, formulae (EJ) H- (EJ) represent one of the most natural recipes 
for the practical construction of the metric. In this sense our present results may be 
formulated, briefly, as a successful application of this recipe to the sequence of matrix 
models as sampled, at even N, by Eq. ffTTj) in Appendix B. 

For our present purposes we decided to study just a one-parametric subset of this 
family. Having skipped the N = 2 model as trivial, let us pick up, for illustration, the 
next, N = 4 example [19] . Once we reparametrize the Hamiltonian H^ b ~. in terms of a 
distance A from the exceptional-point value A^ p ^ = 0, 



-3 VzVT^ 
-VzVT^J -1 2 V / T^A 

v^v/r^A 








-2V1^A 1 

o -V3y/T^x 



(9) 



we reveal that the new parametrization simplifies the spectrum, 



En 



-3y/X, Et = -V\, 



V\, E 3 = 3V\. 



It remains real at any A G (0, oo) = x>( physical \ The comparison of this result with its 
trivial N = 2 predecessor [17] opened in fact the way towards its extrapolation to all N, 
even or odd, as given, in Appendix B, by Eq. (|T2|) . 
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3.1 The simplest metrics and observables: N 



Any two-by-two real-matrix candidate could play the role of a toy model for which 
we would be able to specify the subdomain j)^P hys%cal ' of free real parameters leaving the 
spectrum real. Equally easily we would even construct the related metrics and de- 
termine the Hilbert space(s) in which the matrix would become, constructively, 
tractable as a self-adjoint Hamiltonian [T7] . 

The weak point of such an approach would lie in a vast ambiguity of the transition 
to any larger N by N matrix H^ N \ That's why our present selection of the class of 
Hamiltonians started at the large N. Once we specified the general structure of H^ N ' we 
were able to study its consequences in more detail. 

Let us now return to the first nontrivial dimension N = 2 at which our toy-model 
Hamiltonian reads 

-1 VT^X 

- vt^x i 



H^(X) 



This is precisely the matrix which also appeared, in a slightly different notation and in 
an entirely different physical context, in our preceding paper [25J. Still, the mathematics 
remains the same, including the recommended replacement of A by r = r(A) = y/X. 
This has been shown to facilitate the solution of the conjugate Schrodinger Eq. (jTJ), i.e., 
equivalently, via Eq. (jSJ), the construction of the necessary Hermitizing metric in closed 
form. Secondly, the introduction of the further two apparently redundant abbreviations 
y/l — r := u = u(X) and \/l + r := v = v(X) enabled us to write down the resulting ketkets 
= {u, — v} T and \ipo)) — {~ v,u} T in a particularly compact and symmetric manner 
(notice that the superscripted T means transposition while the ordering of the ketkets is, 
for some inessential reasons, reversed). 

Once we now return to the N = 2 case (and to the construction of the metric) , another 
abbreviation appears also useful, z — uv 
the matrix of metric most compact, 



\J\ — A. This in fact makes the sum (IS|) for 



= 6 (2) (A) 



1 -VT^X 
-VT^X 1 

One might also notice that this matrix becomes diagonal in the far-from-QC limit A — > 1. 





1 -z 




-z 1 
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In opposite direction, any other observable quantity must be guaranteed to be crypto- 
Hermitian, i.e., compatible with the Dieudonne's [7] constraint G^Q = QG. Thus, our 
knowledge of the metric is vital. It implies that any N = 2 candidate for an operator of 
an observable represented, say, by the ansatz 



G 



a b 
c d 



must be restricted by the Dieudonne's constraint yielding the single rule c — b = z(a — d). 
Optionally, we may make the new observable G independent of A (i.e., of z). For this 
purpose it suffices to satisfy the two separate constraints a = d and c = b. In the latter 
case, the ultimate condition of the reality of the eigenvalues of G is trivially satisfied for 
any real d and b. 



3.2 The metric at N = 3 

At the first sight the use of the doublet of ad hoc functions u = u(X) and v = v(X) of A 
seems artificial. Recently we imagined that their appeal need not necessarily be restricted 
to iV = 2 and rewrote the triple-level Hamiltonian in these variables, 

-2 V2uv 
—y/2uv y/2uv 
-V2uv 2 

We felt that this introduces a new symmetry in our present problem. In a backward 
perspective, this idea truly will play the key role in our forthcoming considerations. 

On the level of energies, the idea looks trivial. Indeed, after transition from N = 2 
to N = 3 we merely obtain a robust, coupling-independent new level E — 0. From the 
qualitative point of view, no new insight in the QC mechanism of spectral complexifications 
is gained. The situation only becomes more interesting when we recall formula (jHJ) and try 
to define the metric. This must be preceded by the construction of the left eigenketkets of 
the Hamiltonian H^ 3 '(u,v). After some elementary calculations one obtains the formula 
which can also be found, mutatis mutandis, in Ref. [25], 

M) = ~V2uv, v 2 } T , |Vx» = {~V2uv, 2, -V2uv} T , 
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|Vb)) = {v 2 ,-V2uv,u 2 } T ■ 
The comparison with its A/" = 2 predecessor does not offer any hint for extrapolation to 
N > 3. In a search for some new symmetries, the latter result must be further amended. 

By the trial and error technique we succeeded in verifying that the key to a successful 
amendment should be sought in the homogenization of the individual ketkets, treated 
exclusively as functions of u and v. At iV = 3 the only item which violates such an overall 
principle is the second element of which is equal to 2. Thus, in place of this digit 
we will insert the expression u 2 + v 2 which is identically equal to 2. This seems to be the 
trick. After we insert the amended vectors in Eq. (JS}, we obtain the metric in its highly 
indicative form of the matrix sum 








V2 









e< 3 > = i — z 


y/2 





y/2 


+ z 2 J, 


z = uv 







V2 










where the symbol J denotes the matrix with units along the second diagonal. This matrix- 
polynomial version of the metric at A/" = 3 seems to exhibit all of the expected symmetries. 
Vice versa, the presence of these "hidden" symmetries may be expected to simplify also 
the matrix structure of the metric at the higher dimensions N. 

3.3 Crypto-Hermitian observables at N = 3 

Every unique specification of the metric determines the complete family of the eligible 
observables G. The only constraint is that these matrices must be crypto-Hermitian, i.e., 
compatible with the Dieudonne's crypto-Hermiticity condition G^Q = QG. From the 
ansatz 





i 


-V2z 


z 2 




a 


b 


c 


e = 


-y/2z 


1 + z 2 


-V2z 


G = 


d 


f 
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z 2 


-V2z 


1 




h 


k 


m 



this enables us to deduce the triplet of explicit crypto-Hermiticity conditions imposed upon 
the general candidate matrix G, 

z 2 a — \/2zd + h — c + \/2zg — z 2 m = 0, 
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-V2za + d + dz 2 - V2zh -b + V2zf - z 2 k = 
z 2 b - V2zf + k + V2zc - g- gz 2 + V2zm = . 
The general solution may be written in the following form, 

b = -V2za + d + dz 2 - V2zh + \p2zj - z 2 k, 

c = z 2 a- V2zd + h + z 2 m + 2z 2 h-2 z 2 f - V2z 3 d + V2z 3 k + V2zk, 

g = V2zm + V2zh - V2zf - dz 2 + z 2 k + k. 

Whenever needed, our observables may even be required z-independent (i.e., indepen- 
dent of the A— controlled changes of the Hamiltonian) . The explicit form of the latter 
observables is, naturally, less flexible. Still, they may be shown to possess the following 
three-parametric general form 





a 


d 


h 


F = 


d 


a + h 


d 




h 


d 


a 



The three related eigenvalues 

|a — h, a + h + V2d, a + h — v^2c?| 

should be required real of course. Thus, it is important to conclude that the latter subfamily 
of specific observables F remains conserved and uninfluencedby the occurrence of the global 
quantum catastrophe in the QC limit z — > 1. 

4 The extrapolation pattern 

The technical core of our present message lies in the description of merits of a rather 
counterintuitive replacement of the single parameter A by an apparently redundant pair 
of its two functions u = u(\) and v = v(X). This trick enables us to make use of the 
symmetries hidden in the problem. As a consequence, an efficient extrapolation strategy 
may be developed, simplifying decisively the use of the spectral formula ([8]) for metric via 
the decisive clarification of the structure of the necessary set of the left eigenvectors of our 
conjugate Hamiltonians H^ N \\) at any matrix dimension N. 
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4.1 The first nontrivial case — the confluence of the two com- 
plexifications at N = 2J with J = 2 

At N = 4, our toy model experiences its QC collapse at the boundary of the accept- 
ability interval, i.e., at z — z^ P ^ = ±1 where the quadruplet of the energies completely 
degenerates, E^ P ^ = 0. The parallel QC degeneracy involves also the wave functions. 
They collapse into the same eigenvector in the QC limit as well. The same degeneracy rule 
applies, finally, also to the left eigenvectors of H, i.e., the right eigenvectors of W . 

The latter statement is not so easily checked at N = 4. Indeed, the direct calculations 
using symbolic manipulations (e.g., in MAPLE) offer just the quadruplet of the eigenvectors 



of 



it 



#(4) 

(z) in a rather complicated closed form, 

( 3 VT^z 2 + 9z 2 -12 y/l^z 2 - 12 N 
3 z (-z 2 + 2 + 2 VT^z 2 ) V3 



|V>3»2 
|^3»3 
V l^3»4 J 



\ 



' 1^2 






)>2 


1^2 


»s 


\ 1^2 


)>4 , 


f M 






2 




3 




4 / 



-3z 2 V3(VT^+l) 
3z 3 

( -3z 2 (VT^+l) \ 
zV3 (2 y/T^z 2 + z 2 + 2) 
-2 2 v / 3(v / r 3 ^ + 3) 



/ 



3z 3 



/ 



3z 2 (vT 
rv^ (2^ 



/ 

1) ^ 

z 2 -2) 



! \/3(v / r^ 2 -3) 



V 



3z 3 



/ |-0o»i \ / -3 v 7 ! 3 ^ 2 + 9 z 2 + 12 v 7 !^ 2 



l^o)) 2 

V l^0»4 J 



12 



\ 



-3 z (z 2 - 2 + 2 v 7 ! 3 ^ 2 ) 
3^ 2 v / 3(v / T 3 ^ 2 - 1) 

3^ 3 

The apparently complicated structure of this result is deceptive, and a perceivable com- 
pactification of these formulae is possible. In Ref. [25] we only managed to obtain a partial 



V 
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answer, viz., a solution in a hard-to-extrapolate form 

|V> 3 )) = K, -v^m^, v^m; 2 , -w 3 } T , 
M = {v 7 ^, -(3 + r)«, (3 - r)v, -v / 3mw 2 } t 
A more extrapolation-friendly structure of these formulae is needed. The trial-and-error 
method led us to the success when we applied the identities 3 + r — > u 2 + 2 v 2 and 3 — r — > 
2u 2 + f 2 . This finally produced the extrapolation-friendly result. On this basis we may 
now conjecture the following general formula 



( \1>n-i)) ^ 

|-0AT-2» 



A? 



(10) 



V |-0o» J 

The individual matrix coefficients are assumed diagonal (M {N) (1) = I), bidiagonal (M {N) {2)), 
tridiagonal and "rhomboidal" (^^(3)) etc, ending up with the same antidiagonal Ai^ N \N) - 
J = \fl as above. 

At N = 4, the validity of this conjecture may be confirmed by the brute-force solution 
of Eq. (|7|). The explicit version of the formula for the left eigenvectors is obtained, 

10 
10 
10 
1 
V3 



u 



+ 








2 0^ 



VS 2 








uv 



V3 
s/?> 2 
2 ^ 

\/3 

1 

10 
10 

10 



M 2 V+ 



This result exhibits the expected symmetries as well as the extrapolation-friendly sparse- 
matrix structure of the individual expansion matrices. At any dimension N the knowledge 
of this pattern will decisively facilitate the concrete determination of the numerical values 
of the matrix elements as well as the ultimate use of Eq. ([8]). 
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4.2 A confirmation of the pattern: N 



At iV = 5 our Hamiltonian matrix still fits in the printed page, especially if we abbreviate 

z = uv, 



tf(5), 



u,v) 



-4 


2z 











-2z 


-2 


\/Gz 











— y/Ez 





\/6z 











—a/6 z 


2 


2z 











-2z 


4 



After all of the symbolic manipulations needed we get the expected extrapolations of the 
coefficient-matrices in ( fTOi) at N = 5; accompanied by the "missing", not yet predicted 
two items 



M {5) (2) 



2 

2 Vg 

v/6 Vq 

0^02 

2 



a/6 

3 3 

a/6 4 v/6 

3 3 

Vg 



This observation demonstrates that the move to any higher dimension N becomes easily 
implemented via the sparse-matrix ansatz f fTUj) . In the same spirit, also the construction 
of the metric becomes just a matter of a very routine linear algebra. The only obstacle 
emerges due to the growth of the size of the resulting matrices. As long as they do not not 
fit in the printed page anymore, their elements must be displayed, whenever needed, in a 
suitably compressed form (in a different context, interested readers may find a sample of 
such a compression in Ref. 



5 Summary 

In contrast to the classical Thorn's theory of catastrophes [26], it seems rather difficult to 
formulate the very purpose of its sufficiently satisfactory quantum counterpart, not even 

13 



speaking about its mathematics itself. In this context we described here just one of many 
possible approaches to the problem. 

During the preparatory and purely formal considerations we imagined that from the 
pragmatic, phenomenological point of view, a smooth change of a suitable parameter may 
lead, in many models, to an abrupt complexification of some energy level or levels, i.e., to 
an abrupt loss of their observability status. This is what we decided to call here a quantum 
catastrophe. 

In our text we simulated the QC process (during which the parameter crosses its 
exceptional-point value A^ EP - ) G <92}( p/ls/Sica/ )) via suitable toy models. Our motivation was 
obvious: as long as the textbooks on quantum theory just rarely cover the QC phenomena 
in a systematic manner, the field may be perceived as open to new theoretical develop- 
ments. In parallel, the recent reformulations of the representation theory characterized by 
the use of a nontrivial inner product in the physical Hilbert space appeared suitable for 
the purpose. Last but not least, we felt encouraged by the recent growth of activity in 
experimental physics where the question of relevance of the Kato's exceptional points in 
quantum phenomenology and measurements has been revitalized in several directions [10] . 

An overall mathematical difficulty of the problem (and, in particular, of its more so- 
phisticated, fine-tuned N = 2 J > 4 versions) exposed us to the necessity of choice between 
a phenomenological numerical study of some realistic models (this is the way we choose 
in our recent papers [25]) and an instructive non- numerical description of some carefully 
selected toy models. 

In the present paper we opted for the latter. During our lucky choice of the family of 
models we felt attracted by the unexpectedly friendly nature of their spectra and of their 
geometry in the space of the dynamics-determining parameters (this is an older result 
recollected in Appendix B). In the present text we complemented these observations by 
the discovery of an equally unexpected friendliness of these models from the point of view 
of the systematic construction and extrapolation of their metrics to all dimensions. 

In summary, we believe that our toy models will offer a useful guidance for contin- 
ued research in the field of quantum catastrophes. Not only due to the feasibility of our 
present constructions but also due to the transparency of the matrix structure of their 
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metrics 0^. Indeed, as long as these matrices specify the inner products in the physical 
Hilbert spaces of states 7i^ s \ their compact form opens the way not only towards a facili- 
tated physical probabilistic interpretation of the quantum systems in question but also, as 
we demonstrated, to an unexpectedly transparent matrix structure of the other, generic 
crypto-Hermitian observables. 
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Appendix A. Crypto- Her mitian Hamiltonians 



In Ref. [I] we introduced the notation in which the same state tp of a quantum system in 
question is represented by three alternative ket- vector elements \if)V') of the respective dif- 
ferent Hilbert spaces Hy' with superscripts j = P, F, S. The meaning of these superscripts 
|Jj is given by the following scheme, 



primitive physical space 



simplification ^/ 



unitary equivalence 



F — false space 



hermit ization 



standard physical space 



The role of the arrows is the following. Firstly, the "simplification" arrow means that 
the presumably complicated ket is redefined as the so called Dyson's map of a 

simpler ket \ip^), viz., \tp^) = |^^)- Secondly, the bidirectional "unitary equivalence" 
relationship requirement implies that for non-unitary f2s the Hilbert space must be 
endowed with a nontrivial metric O = Q^Q [2J. Thirdly, the "hermit ization" arrow should 
be read as a replacement of the conventional Hermitian conjugation H — > H' using trivial 
0^) = I by the crypto-Hermitian conjugation H — > H$ := _1 if using unconventional, 
nontrivial = 0^ ^ I. 

One tacitly assumes that the given Hamiltonian H = H(A) has a real spectrum (i.e., 
that A G z>( physi,ca ^) and that the selected metric is a bounded, invertible and positive 
definite operator. In such a dynamical regime there is no real reason for calling our 
Hamiltonian (defined as acting in both of the spaces T-L^'^) non- Hermitian. It is more 
natural to declare a false or manifestly unphysical space. Unfortunately, this space 
is precisely the space in which we make all calculations. Often, it is chosen in the most 
common form L 2 (M) in which the kets are represented by the quadratically integrable 
functions ip(x) where the real variable x does not represent the observable position [30J. 

The latter conventions often become the source of misunderstandings. For this reason, 
the Hamiltonian H (which is, in full compatibility with the first principles of quantum 
theory, Hermitian in its proper and manifestly physical Hilbert space U {s) ) should better 
be called crypto-Hermitian (this emphasizes the not too frequently encountered fact that 
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the correct physical metric is chosen nontrivial, O^ ^ I). 

In the literature, the crypto-Hermitian operators H are also known as quasi-Hermitian. 
In this case one should have in mind the newer definition used in Ref. [2] and not the older, 
more abstract one which has been introduced, by Dieudonne [7j, in the context of pure 
mathematics. 

The majority of physicists who write about the subject think that one should put the 
main emphasis upon the practical aspects of the quantum model in question. One of these 
aspects is that the concrete physical predictions (e.g., the spectra of bound-state energies) 
may still be based on the calculations performed in the friendly-space representation of ip 
(and then just transferred to physical for interpretation purposes). This explains why 
one finds the crypto-Hermitian (i.e., in their proper space HS S \ Hermitian) operators H(X) 
with real spectra still called, in the large number of truly serious and influential papers, 
"non-Hermitian" . 

Naturally, the situation gets changed when the parameter A leaves the physical domain 

V (physical) and when at leagt 

some of the energies become complex. Then, the name "non- 
Hermitian H(\y becomes fully deserved. Indeed, in our three-Hilbert-space scheme of 
Ref. [I], two out of the three Hilbert spaces (viz., / HS S > and H^) simply cease to exist. One 
must find another, open- or sub-system [5] or resonance-theory [B] physical interpretation 
of the quantum system in question. 

At A = \( EP ^ E gj)(phy steal) ^ • e ^ exceptional, quantum-catastrophe values of the 

parameters there also does not exist any reasonable physical interpretation of the physical 
system under consideration. In particular, the metric ceases to exist so that the system 
does not possess any standard quantum interpretation. The measurability status of the 
energies survives (they are still all real) but some vital dynamical information is missing. 
The mechanism of the quantum catastrophe is unspecified. If asked for, it must be added 
via an appropriate enrichment of the model. Typically, such enrichments are in active use 
in magnetohydrodynamics [14J or in laser-physics |31j . 
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Appendix B. Solvable N by N models 



The study of properties of the general N by iV Hamiltonians H^ N ' remains non-numerical 
up to the dimension N = 4, i.e., up to the secular polynomials of the fourth order (cf. |19j). 
At the higher N the spectra are usually studied by numerical or perturbation methods. In 
the latter context the most popular models are the so called anharmonic oscillators, 

1 ... 
3 ... 
5 '•• 



0{g) 



which were chosen, in Ref. [22], as a starting point of a simplification intended to lead to 
non-numerically tractable toy models. 

We restricted our attention to the tridiagonal and antisymmetric perturbations. We 
also truncated our Hamiltonians to iV by iV matrices. With even N = 2J this yielded the 
sequence 



H 



(2) 
(«) 



1 a 
-a 3 



H 



(4) 

(a,b,c) 



1 b 

-b 3 
-a 




a 

5 c 



-c 7 

Finally, we shifted the energy scale and imposed an additional symmetry on perturbations. 
For the resulting set of toy-model Hamiltonians 

-3 b 

— 1 a 



-a 1 



11 



-b -1 a 
-a 1 b 
0-6 3 

the secular equations simplified so that the bound of feasibility of constructive consid- 
erations grew up to iV = 11 [23]. Still, the comparatively large number J = [N/2] of 
variable matrix elements kept the model sufficiently flexible and well adapted to many 
phenomenological needs [T7J [HI [20] . 

In the language of mathematics the main results concerned the conditions of the reality 
of the energy roots. The proofs have been rendered possible by the computer-assisted 
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symbolic manipulations. A brief summary of some technical aspects of these manipulations 
may be found in Ref. [32]. In the second stage of developments we reanalyzed the secular 
equations and sought for the strong-coupling extremes of the J = [N/2]— dimensional 
real domain x>^ phys%cal \ As a result we obtained the series of matrices representing the 
degenerate Hamiltonians at the QC instant, 



H, 



(2) 



QC 



1 1 
- 1 -1 



H, 



(4) 
QC 



3 


Vs 








Vs 


i 


2 








-2 


-1 


V3 








-V3 


-3 



Elements a^ QC \ b^ QC \ . . . were determined, within the Grobner-basis elimination method, 
as roots of a polynomial. Although the degree of this polynomial was quickly growing with 
J (e.g., it was already 17 at J = 4 [32]), we still managed to find and prove the general 
extrapolation pattern. The closed-formula version . was found for all N, 

even and odd (cf. Ref. [22J). 

Although the physical meaning of the model is, naturally, completely lost in its fully 
degenerate QC limit, we specified our next task as the reconstruction of the standard 
probabilistic interpretation of the system in the close vicinity of this singularity. We 
revealed that for such a purpose it makes sense to reparametrize our original family of 
Hamiltonians in terms of a distance A > of the Hamiltonian from its exceptional-point 
extreme at X^ EF "> = 0. At iV = 2 this was easy. We merely replaced the above-mentioned 

(2) 

one-parametric matrix H)J. by its reparametrized alternative 



#S(A) 



1 y/l-A\ 
VT^AJ -1 



depending just on the product A A. At iV = 4 one easily arrives at the less trivial Hamil- 



tonian Hrf B] (X) 



[A,B]\ 



3 v^Vl - A - B\ 2 

V3V1-A-5A 2 1 

-2V1 - A - AX 2 
















2 VI - A - AX 2 
-1 

V3V1- X-BX 2 
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etc. 

In Ref. [23] we managed to prolong the series of these reparametrizations to all N = 2 J 
or N = 2 J + 1. We revealed that in the new parametrization of the vicinity of the QC 
extreme, the geometry of the interior of the domain j)(P h y stcal ) becomes trivial, viz, flat 
and layer-shaped. This means that in QC regime, our J— dimensional domain p(p% slca 
becomes characterized by the mere single inequality at any J. Thus, excluding the first, 
slightly anomalous N = 2 case we obtained the sequence of inequalities 

-/4 < 2A/2 — B < +v\ , N — 4, 

-/il < 6A/2 - AB + C < W & , N = 6, 

-t4 < 20A/2 - 155 + 6C - D < +ul , N = 8, 

etc, with /i 4 = 1/2 and z/ 4 = 2/3, etc. It is worth noticing that the coefficients in these 

/ N-2 \ 

layer-specifying inequalities are just the combinatorial numbers I 

\ m I 

In our very recent application-oriented paper [25J we decided to choose just the sim- 
plest, one-parametric subset of the J-parametric Hamiltonian-matrix series H^' B ,(A) 
of Ref. [23] ■ For the sake of simplicity we selected A = 1 at J = [N/2] = 1 and 
A = B = ... = at J > 2. Under these assumptions, we were able to complete the 
task and to construct the metrics G, in a "brute force" manner, up to iV = 8. 

A strong motivation for our present return to the underlying mathematics may be seen 
in the fact that the corresponding energy eigenvalues were found to form the equidistant 
set at any positive iV and A, 

E n = (2n + l- N)^X, n = 0,l,...,N -1. (12) 

Such an unexpected and important merit of the model appeared to be in a sharp contrast 
with the feasibility limitations to ^ 8 as encountered in Ref. [23]. In this sense, our 
present paper just solves this puzzle and outlines the pattern of extension of the construc- 
tion of the metric to any dimension N. 
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